Polyphyletic Alyssum cuneifolium (Brassicaceae) revisited: Morphological and genome size differentiation of recently recognized allopatric taxa

Alyssum cuneifolium has been recognized as a perennial alpine species growing in five isolated European mountain ranges: the Pyrenees, Western Alps, Apennines, Pirin Mts and Mt Smolikas. Recent molecular systematic studies revealed that the disjunct populations from distant mountains are not closely related and belong to five independent species: A. cacuminum (Spain, Pyrenees), A. cuneifolium (Italy, Apennines), A. flexicaule (France, Western Alps), A. pirinicum (Bulgaria, Pirin Mts), and A. spruneri (Greece, Mt Smolikas). The present study brings the thorough morphometric analysis of the segregated taxa. We found minor morphological differences between them. Whereas A. pirinicum can be clearly distinguished, the other taxa are recognizable only at the level of population means of investigated characters. The morphological similarity of these distantly related species is obviously the result of adaptation to similar high‐alpine scree habitats. It is not clear, however, whether this adaptation is environmentally controlled or whether it is also genetically fixed and whether it reflects parallel evolution towards similar morphotypes. The observed morphological patterns and their assumed correlation with environmental factors are discussed using examples from other Alyssum taxa. Three different ploidy levels have been reported for the species under study. In the present article, we examine variation in relative nuclear genome size. The Alpine and Pyrenean species have larger relative monoploid genome sizes than the Apennine and Balkan ones, probably reflecting the evolutionary history of the group. A nomenclatural account of the study species is presented, and lectotypes of A. cuneifolium and of two other names are selected.     

Chloroplast phylogenomic data support Eocene amphi‐Pacific early radiation for the Asian Palmate core Araliaceae

Traditional phylogenies based on analysis of multiple genes have failed to obtain a well‐resolved evolutionary history for the backbone of the Asian Palmate group of Araliaceae, the largest clade of the family. In this study, we applied the genome skimming approach of next‐generation sequencing to address whether the lack of resolution at the base of the Asian Palmate tree is due to molecular sampling error or the footprint of ancient radiation. Twenty‐nine complete plastid genomes of Araliaceae (17 newly sequenced) were analyzed (RAxML, Beast, Lagrange, and BioGeoBears) to provide the first phylogenomic reconstruction of the group (95% of genera included). As a result, the early divergences of the Asian Palmate group have been clarified but the backbone of its core is not totally resolved, with short internal branches pointing to an ancient radiation scenario. East Asia is inferred as the most likely ancestral area for the Asian Palmate group (from late Paleocene to Eocene) from which early colonization of the Neotropics is inferred during the Eocene. The radiation of the core Palmate group took place during the late Eocene, most likely in the context of the Boreotropical hypothesis. Recurrent episodes of southward migration (to the tropics) coupled with northern latitude local extinctions (promoting geographic isolation of lineages) followed by northward expansion (promoting contact of lineages that erased the trace of preceding geographic isolation) are hypothesized to have linked to the alternation of the cold and warm periods of the Eocene.     

Polyploidy does not control all: Lineage‐specific average chromosome length constrains genome size evolution in ferns

Recent studies investigating the evolution of genome size diversity in ferns have shown that they have a distinctive genome profile compared with other land plants. Ferns are typically characterized by possessing medium‐sized genomes, although a few lineages have evolved very large genomes. Ferns are different from other vascular plant lineages as they are the only group to show evidence for a correlation between genome size and chromosome number. In this study, we aim to explore whether the evolution of fern genome sizes is not only shaped by chromosome number changes arising from polyploidy but also by constraints on the average amount of DNA per chromosome. We selected the genus Asplenium L. as a model genus to study the question because of the unique combination of a highly conserved base chromosome number and a high frequency of polyploidy. New genome size data for Asplenium taxa were combined with existing data and analyzed within a phylogenetic framework. Genome size varied substantially between diploid species, resulting in overlapping genome sizes among diploid and tetraploid spleenworts. The observed additive pattern indicates the absence of genome downsizing following polyploidy. The genome size of diploids varied non‐randomly and we found evidence for clade‐specific trends towards larger or smaller genomes. The 578‐fold range of fern genome sizes have arisen not only from repeated cycles of polyploidy but also through clade‐specific constraints governing accumulation and/or elimination of DNA.     

海南产桶形芋螺线粒体基因组DNA提取条件的优化

提取海南产桶形芋螺线粒体基因组完整DNA (mtDNA),并对提取条件进行优化。以桶形芋螺腹足肌肉、毒腺和肝胰脏三个不同组织为材料,分别采用改进高盐沉淀法、细胞器/磁珠法和试剂盒提取三种方法,提取桶形芋螺mtDNA,并利用琼脂糖凝胶电泳和紫外分光光度计对提取mtDNA的纯度和浓度进行测定。以coxⅠ-rRNA小亚基基因和α-芋螺毒素基因设计引物,通过PCR反应来确证所提取的DNA确实是mtDNA。试剂盒法提取肝胰脏、高盐沉淀法提取肝胰脏和腹足肌肉组织这三种方法的产率很高,分别为44.4μg/mg、43.3μg/mg和32.6μg/mg。A260/280比值表明,改进高盐沉淀法提取毒腺和腹足肌肉组织,细胞器磁珠法提取腹足肌肉组织的mtDNA纯度很高。综合比较,采用改进高盐沉淀法,利用桶形芋螺腹足肌肉组织所提取的mtDNA产率高、质量好、纯度高。高质量芋螺mtDNA的获取为利用分子生物学方法对芋螺进行遗传进化分析和系统分类提供了基础。     

Finding stories in noise: Mitochondrial portraits from RAD data

Mitochondrial DNA (mtDNA) has formed the backbone of phylogeographic research for many years; however, recent trends focus on genome‐wide analyses. One method proposed for calibrating inferences from noisy next‐generation data, such as RAD sequencing, is to compare these results with analyses of mitochondrial sequences. Most researchers using this approach appear to be unaware that many single nucleotide polymorphisms (SNPs) identified from genome‐wide sequence data are themselves mitochondrial, or assume that these are too few to bias analyses. Here, we demonstrate two methods for mining mitochondrial markers using RAD sequence data from three South African species of yellowfish, Labeobarbus. First, we use a rigorous SNP discovery pipeline using the program stacks , to identify variant sites in mtDNA, which we then combine into haplotypes. Second, we directly map sequence reads against a mitochondrial genome reference. This method allowed us to reconstruct up to 98% of the Labeobarbus mitogenome. We validated these mitogenome reconstructions through blast database searches and by comparison with cytochrome b gene sequences obtained through Sanger sequencing. Finally, we investigate the organismal consequences of these data including ancient genetic exchange and a recent translocation among populations of L. natalensis, as well as interspecific hybridization between L. aeneus and L. kimberleyensis.     

An evaluation of sequencing coverage and genotyping strategies to assess neutral and adaptive diversity

Whole genome sequences (WGS) greatly increase our ability to precisely infer population genetic parameters, demographic processes, and selection signatures. However, WGS may still be not affordable for a representative number of individuals/populations. In this context, our goal was to assess the efficiency of several SNP genotyping strategies by testing their ability to accurately estimate parameters describing neutral diversity and to detect signatures of selection. We analysed 110 WGS at 12× coverage for four different species, i.e., sheep, goats and their wild counterparts. From these data we generated 946 data sets corresponding to random panels of 1K to 5M variants, commercial SNP chips and exome capture, for sample sizes of five to 48 individuals. We also extracted low‐coverage genome resequencing of 1×, 2× and 5× by randomly subsampling reads from the 12× resequencing data. Globally, 5K to 10K random variants were enough for an accurate estimation of genome diversity. Conversely, commercial panels and exome capture displayed strong ascertainment biases. Besides the characterization of neutral diversity, the detection of the signature of selection and the accurate estimation of linkage disequilibrium (LD) required high‐density panels of at least 1M variants. Finally, genotype likelihoods increased the quality of variant calling from low coverage resequencing but proportions of incorrect genotypes remained substantial, especially for heterozygote sites. Whole genome resequencing coverage of at least 5× appeared to be necessary for accurate assessment of genomic variations. These results have implications for studies seeking to deploy low‐density SNP collections or genome scans across genetically diverse populations/species showing similar genetic characteristics and patterns of LD decay for a wide variety of purposes.     

Chromosome‐level reference genome of X12, a highly virulent race of the soybean cyst nematode Heterodera glycines

Soybean cyst nematode (SCN, Heterodera glycines) is a major pest of soybean that is spreading across major soybean production regions worldwide. Increased SCN virulence has recently been observed in both the United States and China. However, no study has reported a genome assembly for H. glycines at the chromosome scale. Herein, the first chromosome‐level reference genome of X12, an unusual SCN race with high infection ability, is presented. Using whole‐genome shotgun (WGS) sequencing, Pacific Biosciences (PacBio) sequencing, Illumina paired‐end sequencing, 10X Genomics linked reads and high‐throughput chromatin conformation capture (Hi‐C) genome scaffolding techniques, a 141.01‐megabase (Mb) assembled genome was obtained with scaffold and contig N50 sizes of 16.27 Mb and 330.54 kilobases (kb), respectively. The assembly showed high integrity and quality, with over 90% of Illumina reads mapped to the genome. The assembly quality was evaluated using Core Eukaryotic Genes Mapping Approach and Benchmarking Universal Single‐Copy Orthologs. A total of 11,882 genes were predicted using de novo, homolog and RNAseq data generated from eggs, second‐stage juveniles (J2), third‐stage juveniles (J3) and fourth‐stage juveniles (J4) of X12, and 79.0% of homologous sequences were annotated in the genome. These high‐quality X12 genome data will provide valuable resources for research in a broad range of areas, including fundamental nematode biology, SCN–plant interactions and co‐evolution, and also contribute to the development of technology for overall SCN management.     

Phylogenomics using low‐depth whole genome sequencing: A case study with the olive tribe

Species trees have traditionally been inferred from a few selected markers, and genome‐wide investigations remain largely restricted to model organisms or small groups of species for which sampling of fresh material is available, leaving out most of the existing and historical species diversity. The genomes of an increasing number of species, including specimens extracted from natural history collections, are being sequenced at low depth. While these data sets are widely used to analyse organelle genomes, the nuclear fraction is generally ignored. Here we evaluate different reference‐based methods to infer phylogenies of large taxonomic groups from such data sets. Using the example of the Oleeae tribe, a worldwide‐distributed group, we build phylogenies based on single nucleotide polymorphisms (SNPs) obtained using two reference genomes (the olive and ash trees). The inferred phylogenies are overall congruent, yet present differences that might reflect the effect of distance to the reference on the amount of missing data. To limit this issue, genome complexity was reduced by using pairs of orthologous coding sequences as the reference, thus allowing us to combine SNPs obtained using two distinct references. Concatenated and coalescence trees based on these combined SNPs suggest events of incomplete lineage sorting and/or hybridization during the diversification of this large phylogenetic group. Our results show that genome‐wide phylogenetic trees can be inferred from low‐depth sequence data sets for eukaryote groups with complex genomes, and histories of reticulate evolution. This opens new avenues for large‐scale phylogenomics and biogeographical analyses covering both the extant and the historical diversity stored in museum collections.